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1 Introduction 

In the last two decades, the Econophysics community has 
deeply investigated the empirical features of historical time 
series emerging in different financial contexts, ranging from 
high to low frequency data from the stock exchange to the 
markets in which bonds, FX rates, commodities, energy, 
futures, options and many other instruments are traded. 
Numerous empirical evidences have emerged, sometimes 
confirming results from the financial-econometrics com- 
munity, or finding new stylized facts and opening new 
fields of research [l][2l|3l|4] . A well known example is given 
by the observation that the intra-day or daily stochastic 
dynamics of stock prices significantly deviates from the 
behaviour predicted by the geometric Brownian motion. 
This finding dates back to the work of Mandelbrot 0, 
whose attention was mainly focused in recognizing realiza- 
tions of stable processes, and to the analysis of Fama [B] 
concerning the long tailed nature of the Dow Jones In- 
dustrial Average single components. The statistics of re- 
turns has been recently rediscussed in different flavours 
and its modelling has considerably grown. Very hetero- 
geneous models, able to reproduce the degree of asym- 
metry and the excess of kurtosis of the measured distri- 
butions, have been proposed. A non exhaustive list in- 
cludes approaches developing from specific distributional 
assumptions, as it is the case of the Levy flights [l5,7,,8^, 
the Generalized Student-t or Tsallis distributions |2l[9llT0] 
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and the exponential one and it can be extended to 
more sophisticated models capturing the stochastic nature 
of the volatility, see [T2l[T3l[T4l[T5] and references therein, 
the power law scaling of its autocorrelation function and 
the leverage effect pifTB] and the multi-fractal properties 
of historical time series [TTIITS]. The level of complexity 
grows increasing the dimensionality of the problem. The 
leptokurtic nature of returns distributions has been shown 
also for financial indexes, which are linear combinations 
of plain stocks, such as the Standard & Poor's 500 [TJ 
[H] or the NIKKEI (ID]. From the point of view of the 
Central Limit Theorem, this result was not expected: be- 
cause of the great number of components, the convergence 
towards the Gaussian regime would be ensured also for 
heavy tailed stock returns with finite variance. Moving 
back from the index to the components has proved that 
the single assets themselves share the same power law scal- 
ing, with tail index close to 3 [H]. Indeed, this fact would 
be responsible for an ultra slow convergence to Normality 
for independent stocks; however the existence of strong 
correlations among assets is a well known evidence which 
violates the assumptions of the Central Limit Theorem. 

In this paper we investigate a completely different mech- 
anism at work, even when the single factors governing 
the aggregate behaviour are assumed to be Gaussian. We 
move from the analysis of financial indexes to the case 
of portfolios where non linear instruments, such as op- 
tion contracts, induce deviations from Normality, which 
are no more a consequence of the microscopic dynamics. 
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In the following sections we will review the standard Delta 
Gamma Normal (DGN) approach to the problem of risk 
management for non linear portfolios and we will discuss 
a novel analytical methodology and its efficient numerical 
implementation. In particular, our aim is to evaluate two 
standard risk measures, the Value-at-Risk (VaR) and the 
Expected Shortfall (ES) in the framework of the DGN ap- 
proximation. For an introduction we refer the reader to 
the standard financial literature [22l[23lf24] . As far as the 
single asset case is concerned, many different approaches 
can be adopted. A widely exploited technique in risk anal- 
ysis, requiring no assumptions on the distribution of the 
data consists in using past returns realizations to compute 
risk exposure levels, which are associated to the empirical 
quantiles of the distribution. Despite its simplicity and 
generality, such approach suffers from the typical draw- 
backs related to finite size effects. On the other hand, 
if an assumption is made on the data distribution, var- 
ious parametric approaches are available, both in a fre- 
quentist and Bayesian framework, often leading to ana- 
lytical or semi-analytical expressions for VaR and ES p51 
. Beyond the standard results for Normally distributed 
asset returns, in Econophysics risk measures have been 
evaluated with heavy tailed probability density functions 
(PDF), see for example [2] and approaches based on gen- 
eralized Student-t [10], Tsallis [26] and Truncated Levy 
random variables [IZ]. Equivalent estimates to those for 
i.i.d. Student-t variables can also be recovered in a differ- 
ent approach exploiting Bayesian Product Partition Mod- 
els for Normal independent but not identically distributed 
returns 1281 . Furthermore, advanced approaches based on 
the so-called downside risk and applied to the specific case 
of hedge funds can be found in [^150] . In this paper, 
we follow an alternative approach based on the general- 
ized Fourier representation of the PDF, already developed 
in [57] for the single asset case. Fourier techniques for VaR 
evaluation have been introduced since the work of Rou- 
vinez |31j , and recently Martin |32| provided an extension 
of this framework, obtaining expressions also for the ES 
and its sensitivities. However, despite being in semi-closed 
form, they are not straightforwardly suited to numerical 
evaluation but require a saddle point approximation. On 
the contrary our analytical formulae can be easily and 
efficiently evaluated by means of standard trapezoidal in- 
tegration or Fast Fourier algorithms. 



This article is organized as follows. In section |2] we re- 
view the portfolio model we consider, outlining the avail- 
able analytical information (its characteristic function and 
the corresponding asymptotic tail behaviour of the PDF) . 
In that section we derive the semi-analytical expressions 
for the VaR, ES, and their sensitivities with respect to the 
model parameters, which represent the original contribu- 
tion of this work. In section [3] we compare the numerical 
results obtained through our expressions with those from 
the Monte Carlo simulation of synthetic portfolios. Even- 
tually, the relevant conclusions are drawn in section [4j 



2 A first step beyond the Normal behaviour 

A realistic financial situation usually involves portfolios 
containing relevant quantities of non linear (options) in- 
struments. A description of these portfolios in terms of a 
linear composition of risk factors is inadequate, especially 
when assessing the market risk exposure. An improvement 
consists in taking a second order expansion of portfolio 
variations in the risk factors. In particular, when these 
are assumed to be Normally distributed, one obtains the 
DGN model. Since this model has become a benchmark 
framework in the common practice of portfolio risk man- 
agement, in the present work we do not discuss the limits 
of applicability of the DGN approximation, which are an- 
alyzed in [lllSSlISSj. 



2.1 The Delta Gamma Normal approach 

According to the DGN model, the portfolio price change V 
over a given time horizon At is described by the quadratic 
form 

V ^0 + A^x + l-x^rx , (1) 



where X is the iV-dimensional vector of the risk factors 
which are responsible for the portfolio fluctuations. Here, 
the variation is defined a.s V — W — Wq, where Wq de- 
notes the value of the portfolio at the present time t — to, 
while W is the corresponding value at t = <o + At. The 
DGN model assumes X to be drawn from a multivariate 
Gaussian distribution A/'(0, U) with zero mean and covari- 
ance matrix S. In equation ([T]) € M and A e are 
constants, and so is the real symmetric N x N matrix F 
which accounts for possible non linearities. 

By solving the generalized eigenvalue problem 



CC' ^ S and C' rC^A 



(2) 



with A — diag(Ai, . . . , Aat), equation ([I]) can be conve- 
niently rewritten as 



V = 



N 

E 



S-Y- -I -Y" 



(3) 



where 5 = A and X — CY. The l^'s now represent 
independent standard Gaussian variables, and the corre- 
lation structure between the actual risk factors contained 
in X is now spread across the new S and A parameters 
through relations (2|. 

Interestingly, the moment generating function of the 
random variables appearing in equation p| has been com- 
puted in [33] as a Gaussian integral, yielaing 



E[e' 



,^^(5,y.-H(A,y.^)/2)i 



: exp 



2 1 - 



.(4) 

w G M, and since this is a holomorph function in the neigh- 
bourhood of 0, the result can be extended to the complex 
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domain. So, the characteristic function / of (|3]) can be 
computed: 

/H=E[e-^]=e'^"n— ^=exp(--^^l . 

(5) 

Risk assessment applications typically require a careful es- 
timation of the tails. Even though not required by the VaR 
and ES evaluation presented in this paper, we devote the 
last part of the present subsection to review the analytical 
results for the tails of the PDF of the DGN model, which 
have been fully characterized by Jaschke et al. in |34| . In- 
deed, let us suppose that the eigenvalues have been 
ordered in increasing order, and let n < A'^ be the number 
of distinct values. Let ik be the highest index of the fc-th 
distinct eigenvalue (so that A^^ < . . . < Ai„) and let mk 
be its multiplicity. Then we can define the following useful 
quantities for k — 1, ■ ■ ■ ,n 

j=ifc_i+i 

The left tail behaviour of the PDF p{V) is determined by 
the sign of the lowest eigenvalue A* = A^^ . When A* < 0, 
the left tail displays an exponentially damped power law 
decay, whose rate is actually given by A*. More precisely, 
when V — > — oo, p{V) has the following asymptotic tail 
behaviour 

p{V) = c((5. A) \vr e^/^*+"iVl2WA-| 

x{l + Oil/V\V\)) , (6) 

where 

r (mi - 3)/4 if ai ^ 
\mi/2-l if ai = . 

When A* = 0, the left tail is characterized by an asymp- 
totically Gaussian scaling: 

piV)=d{6,\) |V^-For5^2-^"^'e-(^-^'')'/(2«?) 

x{l + 0{l/\V-Vo\)) , (7) 

where Vq = 9 — Y^^=2^'k/i'^^ik)- I* worth mentioning 
that, both in this case and the previous one, the contri- 
bution coming from the power law factor may not be neg- 
ligible when estimating relevant quantiles; moreover, in 
this context deviations from Gaussianity arise as a conse- 
quence of the inclusion of second order terms in the port- 
folio variation ([T]), and not of the assumption of a fat tailed 
distribution for the risk factors. 

Lastly, when A_* > the PDF is zero for V < Vinf , with 
Vini — 9 — J2k=i ^i/(2Aifc), and it decays as a power law 
in the hmit V V^'^f. 

piV) = kid, A) [V - Vi^if/^-' (1 + 0{V - Vinf)) . (8) 

In equations ([T]), ([S]), c{S,X) = c{Si=i^,,,^n, K=i,...,n), 
d{S,X) and k(d, X) are constants depending on the set of 
parameters in use. 



As far as risk estimation is concerned, the left tail de- 
cay is what we are interested in; however, previous con- 
siderations apply to the right tail in an antithetic way de- 
pending on the sign of the highest eigenvalue Xi^ . When 
Xi^ > (Ai^ = 0) the right tail is exponential (Gaus- 
sian) in the limit V — >■ +00; otherwise, the support of 
p{V) is limited from the right and for V V^~p it scales 
as a power law with exponent N/2 — 1 and Vsup ~ 9 — 
ELi^^/(2A,J. 

The central moments of the distribution of V can be 
also explicitly derived, and the first four of them read [33l 

mm- 

111 = 9 + ]^tv[rs] 

112 = A^SA + itr[(rr)2] 

/i3 = 'dA^srsA + ii[{rsf] 

/i4 i2A^ sirsfA + 3 iT[{rsf] + iixl . 

From these relations explicit expressions for the skewness 
C — [SJ^'^ and the kurtosis k = /X4//i|— 3 can be derived, 
and one can easily check that whenever r' = we have ^ = 
and K = 0. This is coherent with the fact that equations 
([1]) and ([3]) define a Gaussian portfolio model when their 
quadratic terms are set to zero, i.e. when F — Q. From this 
point of view, the possible presence of asymmetries or non 
Gaussian tails in the PDF of the DGN model stems from 
the non linear terms in ([I]) and ([3]). 

2.2 Formulae for risk estimation 

The risky nature of a portfolio can be accounted for via the 
well-known VaR estimator, which gives the potential loss 
(negative variation) over the time horizon At that could 
be exceeded with probability equal to the significance level 
V* ■ However the VaR is known to suffer from two main 
drawbacks, the lack of subadditivity and of information 
about the average potential loss when the VaR threshold 
is exceeded. Both of them are overcome by introducing an 
alternative and coherent measure [M] known as the ES. 

To obtain semi-analytical expressions for the risk mea- 
sures, the steps followed in [27] are substantially replied. 
We fix a significance level V* € (0, 1) and we define the 
VaR A* as 

V* = I AY p{Y) 

I— A' r-Wo 

= / dV p{V) - dV p{V) . (9) 

In this expression —Wq represents the maximum possible 
loss over At. Compared to [27], this framework involves 
portfolio variations instead of logarithmic returns and this 
leads to the finite lower bound of integration in the first 
line of equation (|9]). Since we are modelling the risk fac- 
tors as Gaussian random variables, p{V) could have un- 
bounded support; nevertheless, under the assumption that 
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the DGN p{V) is a good approximation to the true (un- 
known) PDF, we expect —Wq to be on the far left tail 
of p(V) and the second term of equation ^ to be negli- 
gible [j For this reason, in the following expressions, we 
always imply the limit —Wq — > — oo. 

We now represent p{V) in terms of its generalized 
Fourier transform 



and, recalling the representation (111 oi V* 
pression can be simplified in 



the last ex- 



E*{V*) = A* 



ttV* 



-Re 



+00 



{oj + ij/)2 



(12) 



p(V) = 



1 

2^ 



+ CXD + 1Z-' 



d0 /((/.) e 



-i4,V 



(10) 



where the axis of integration in the complex plane has to 
be chosen parallel to the real axis with imaginary part 
z/ belonging to the strip of regularity (:^_,j^+) of / [371 
[5H] . The boundaries of the strip of regularity are deter- 
mined by the singular points of / closest to the origin. If 
any, singularities must be purely imaginary, and, for the 
present case, a quick analysis shows that they are equal 
to {—\/Xi^, k = 1, . . . ,n}. Using the expression (10) and 



2.3 Sensitivities 

We are not only interested in estimating VaR and ES, but 
also in their sensitivity with respect to the parameters of 
the DGN model. For a portfolio which is a linear compo- 
sition of risk factors. A* and E* are 1-homogeneous func- 
tions [321139] . Indeed, when F — 0, we have A = 5 and, 
posing = for simplicity, the portfolio VaR reduces to 
the well known expression 

A* = -VWlTS V2 crf-i(2P* - 1) , 



switching the integration order, equation ^ becomes where erf ^ is the inverse of the error function. The first 



•p* 



27r 



-l-CJO+ii^ 



-A* 



00-t-i^ 



and the convergence of the innermost integral is guaran- 
teed if we restrict ly S (0,1^"*"). With this choice, the pre- 
vious expression readily reduces to 



•p* 



J ^+00+11^ 
Re 

TT 



+ 00 P ( 



AujA* 



(11) 



where (jj = u + iv. From the previous discussion, we can 
determine the value of v'^ depending on that of A*: 

- when A* < 0, then iy+ = \l/X*\; 

— when A* > 0, then i/^ = +00. 

To evaluate the ES the considerations already exposed for 
VaR evaluation still apply. Starting from its definition and 
using again equation (10), we can write 



E*(V*) - 



1 



2TrV 

^-yA 

TXV* 
-vA' 



-A* 

AV V p{V) 

— 00 

d0 /(0) 

— cxD+if 

-f-00 



iA* 



-Re 



+ 



ttV* 



-A*Re 



(w + 11^)2 



^ The second term in equation ^ can be evaluated numer- 
ically given Wo- Should it be not negligible, the expressions 
presented in this paper are extended in a straightforward way 
by taking into account the surface terms from integrations. 



derivatives of VaR read 

OA* ^ jSS), 

and the homogeneity condition is readily verified 



V2 erf"^(27'* - 1) 



N 



A* 



k=l 



dA* 



It can be shown that an analogous relation holds for the 
E* . The main consequence is that for linear portfolios, the 
knowledge of the first derivatives allows, just by itself, for 
a complete reconstruction of the risk measures. Thus, risk 
measures response with respect to shocks in the Sk weights 
defining the portfolio composition can be fully character- 
ized with no need of computing higher order derivatives. 

On the other hand, in this work we consider the effects 
due to the presence of non linear instruments; in this sce- 
nario, a complete description of the market risk exposure 
would require the knowledge of higher order derivatives. 
Nevertheless, the first order ones still provide crucial infor- 
mation about the sensitivity of risk measures when shock- 
ing the portfolio weights. For this reason, here we report 
the expressions of the first order derivatives only, remind- 
ing the reader that higher order terms can be computed 
going back through the same passages. 

Equation ( 11 ) is now differentiated with respect to the 



generic parameter that, for the sake of tidiness, we address 
as (3 

Since we aim at evaluating the response of the risk mea- 
sures with respect to shocks to the portfolio parameters 
while keeping V* fixed, we obtain 







dV* 



-Re 



+00 



duj ( df{uj + 



ly — icu \ dp 
dA*\ 



+{iuj - iy)f{Lo + iv) 



dp 
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Exploiting the previous equality, and recalling equation ( 10 1, 
the VaR sensitivities read 



dA* 



Re 



Re 
e" 



+00 duj df{ui+iu) iuiA' 



r+00 
Jo i'- 



iui dp 





vA' 



do; /(w + iv) ei" 

dw df{uj + iu) 
V ~\u. 



(13) 



TTp{~A 



-Re 



9/3 



Finally, using again pi|, differentiation of equation (12 1 



gives us the following Fourier representation of the ES 
sensitivities 



dE*{V*) _ dA' J ^ e 



dp 



8(3 [' ttV 

+ 00 



X Re 







■nV* 

^-vA' 

ttV* 



Re 



-Re 



+ 00 



do; ^f{l^J^lv) ^i^^. 



(cj + iz/)2 a,9 

dw a/(a; -f iz^) 



a/3 



(14) 



where the final equality holds because the expression in 
the curly brackets vanishes. Finally we compute the deriva- 
tives of the characteristic function appearing in ( 13 1 and 
(14). Using equation ([s]) we obtain 



5iW 
89 

dm 

85, 

dm 

dXi 



1 - iA,: 
2(1 - iA, 



/(0) , 



1 - 



1 - iA,, 



I{4>) 



(15) 



Incidentally, we notice that all the previous relations, valid 
for the DGN model, are linear in / and thus equations 
and ( 14 ) share a similar structure with equations 
and (n2f . 



3 Numerical results 

In this section we detail the numerical results of our Fourier 
approach, computing A* , E* and the sensitivities for syn- 
thetic portfolios corresponding to the possible cases de- 
pending on the value of the smallest eigenvalue A*. We 
check the semi-analytical estimates with the ones obtained 
by Monte Carlo simulation of the portfolio values. Below, 
a brief summary of the numerical setup for both the semi- 
analytical approach, proposed here, and the standard his- 
torical simulation one is given. 

We do not address the problem of estimating the co- 
variance matrix S from real portfolio data. Nevertheless, 
it is worth noticing that the usual parametric Maximum 
Likelihood estimators of dispersion perform reasonably 
well in the limit of an infinite number of observations for 



the risk factors. When the length of the time series T is 
larger than the number of risk factors iV, it is also possible 
to exploit filtering procedures .42,43j , such as hierarchi- 
cal clustering, to extract the relevant information from 
correlation matrices, retaining only the statistically sig- 
nificant correlations. These techniques, besides reducing 
the dimensionality of the problem, can give crucial infor- 
mation for decision processes such as portfolio allocation. 
Under the thermodynamical limit oi T,N — > +00, with 
T/N fixed, the random matrix theory has also proved 
to be a useful tool to better understand the statistical 
structure of empirical correlation matrices [44(l45] . On the 
other hand, when the sample is small, especially when T 
is smaller than N , which is common for large portfolios, 
the error affecting the Maximum Likelihood estimators is 
large. In this situations, it is possible to resort to better 
performing estimators, such as, for instance, the Ledoit 
and Wolf shrinkage estimator of dispersion j46ll47] which 
is a weighted average of the sample covariance S and a 
target diagonal matrix C 

= {l-a)S + aC , 

where C — (1/rn) X^I'Li Vk^i Vk are the sample estimations 
of the eigenvalues of S and a is the optimal shrinkage 
weight depending on S, C and Xt=i^...,T- 



3.1 Numerical setup 

Fourier inversion. The complex integrals involved by our 
expressions for risk measures and sensitivities are evalu- 
ated by means of trapezoidal integration; this requires to 
split explicitly the real part of the integrand, being an 
even function of oj, in the two terms proportional to the 
sine and cosine functions. For instance, evaluation of V*, 



equation (111, reduces to compute the following two real 



Fourier integrals 



-vA' 



■p* 



+00 



Re 



Im 



cos(a;Z\*) 
sin(ti;Z\*) 



Each one is solved by means of standard adaptive rou- 
tines]^ similar decompositions apply to the other quanti- 
ties of interest. 

Historical simulation. A sample of Tmc realizations of the 
risk factors X ^ N{0, E) is drawn and, correspondingly, 
the time series of the portfolio variations Vt=i^...^TMc 
computed. If Vt are the entries of Vt sorted in ascending 
order, and assuming t* = Tmc x ^* to be integer, the 
historical VaR at the significance level 7^* is defined as|^ 



A*^ = -Vt^ . 



(16) 



^ See the routine QAWFE of the library QUADPACK at 
http:/ /www. nethb. org. 

If Tmc X P* is not integer, A* can be defined as the average 
(— K — Ki)/2, with s, u being the two integers closest to t* . 
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The ES is obtained as the average on the left tail of the 
empirical distribution 



(17) 



Confidence intervals for VaR are obtained from a basic 
result of order statistics. Let us consider the sample of 
Tmc i-i-d. deviates Vi and indicate with Qp the p-th per- 
centile of their distribution; then, the probability of Qp 
being enclosed in (Vt- , Vt+) is given by the following sum 
of binomial probabilities [40] 



k=t- 



P{V,- <Qp< V,.) = {^f)p'{^-pf"^-' • 



(18) 

To find the confidence interval associated to Z\y for a given 
Confidence Level (CL), we find indices satisfying 



P{Vt-- 



< 



<Vt+)<CL< P{Vt- < 



where Q-p* — —A* and the choice between possible dif- 
ferent pairs (t~,t+) satisfying the previous inequality is 
made requiring the confidence interval to be as symmetric 
as possible around A'^. With these positions, Historical 
VaR is estimated by 



(^h) 



+5+ 

-5- 



-Vt- - At 



S- = At 



(19) 



t+ 



with confidence level CL. 

Since E* is a monotonously increasing function of A* , 
a lower and upper bound for ES are easily obtained by 
evaluating the average in equation (17) for t* = t~ and 
t+ respectively. At the same CL as above, we estimate 



t 




(20) 



The sensitivities of Zi^ are evaluated by approximating 
its derivatives with respect to the remapped parameters 
of equation ^ , with the finite difference formula 



{dpA* 



fdA* 



V dl3 



At,iP + Ap) - At,ip ~ AP) 
2AP 



(21) 

where Z\y(/3±Z\/3) denotes the historical estimates for the 
VaR obtained simulating the DGN portfolio after giving 
a positive/negative shock to the parameter j3 while keep- 
ing fixed all the others. The error affecting {dpA*)-^ is 
obtained by linear propagation of the Monte Carlo errors 
on Z\j^(/3±Z\/3) through equation ( plj ). An analogous pro- 
cedure is applied to estimate {dpE*)^ and its confidence 
interval. 



The results discussed in the following sections have 
been obtained simulating Tmc = 10^ synthetic portfolio 
scenarios for = 15 risk factors; for simplicity, we simu- 
lated the remapped factors Y ^ N{0, 1) and the values Vt 
were constructed from equation (|3| after arbitrarily fixing 
6, S and A. Attention has to be paid here to freeze the 
values of Yj, so that re-evaluation of the portfolio with 
shocked parameters is always carried out based on the 
same sample. More details about efficient algorithms to 
generate the scenarios under a quadratic approximation 
of the portfolio losses can be found in [41 . The semi- 
analytical cur ves h ave be en obtain ed b y Fourier inversion 
of equations (111, (12), (13 1 and ( [T4| ), the characteristic 
function and its derivatives as given by ^ and ( 15 ). Given 
a grid of A* values, the corresponding V* , E*^~U^A* and 
dpE* are computed and compared with those from the his- 
torical simulation. If one is interested in A* corresponding 
to a precise spot value of V* ^ inversion of equation (fTTl) 
has to be nested in a root finding procedure. Alternatively, 
the grid of A* can be tuned in such a way to obtain V* 
values close enough to the required one. 



3.2 The lowest eigenvalue is negative 

In order to discuss the case A* < 0, we fix the DGN param- 
eters as follows: 9 = 0, 5i=i ... 15 = 1, A* = A^^i ... 5 = —2, 
A,=6... ,9 = 1 and A,=io,.. ,15 = 2 (CASE 1). 

The correspondence 7'*-VaR and 7^*-ES as obtained 
from equations ( 11 1,( 12 ) is shown in Figure[l] As expected 
from their definitions, E* estimates are always larger than 
Z\*, the shift vanishing in the limit of V* approaching 
0. The sensitivities for this set of parameters are illus- 
trated in Figure [2] and Figure [3j With the exception of 
13 = 9, the derivatives of f{(j)) always depend on the pair 
{Si, Xi); for this reason the only cases we need to consider 
are f3 = 0, (5i,6,iOj Ai^e.io- The sensitivities with respect to 
the central value 9 are identically equal to —1; this re- 



sult can be verified directly from equations (13),(14) after 



substitution of dgf{4)) and reminding the definition of V* 
and p{A*), and the corresponding curves are not shown 
for the sake of clarity. 

All the curves are superimposed to the points obtained 
through historical simulation for the interesting values 
of V* = 0.001, 0.005, 0.01, 0.02, 0.03, 0.04, 0.05, along with 
the corresponding error bars for CL = 98%; central points 
and error bounds correspond to the definitions (19), (20) 
for A* and E* and to their extensions for the sensitivities. 
The latter exhibit larger uncertainty, due to the propaga- 
tion of the statistical error, which is magnified by the finite 
difference formula (21) by a factor of order 0{1/A(3). It 



is also clear how this uncertainty increases for smaller V* , 
due to the decreasing size of the sample on the left tail of 
the distribution of the generated values Vt- On the whole, 
the full agreement with the Monte Carlo outcomes proves 
the effectiveness of our analytical results and the reliabil- 
ity of their numerical implementation. 

As discussed in section [2j in this case the distribution 
function p{V) has both left and right exponential tails. 
This scaling is confirmed by Figure |4| where the left tail 
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0.1 0.5 1 2 3 4 5 

Fig. 1. (Color online) CASE 1 (A* < 0): VaR and ES. The 

semi-analytical curves are compared to the points obtained via 
historical simulation for the values of V* usually considered in 
practice. 
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4^ 
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3 
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Fig. 2. (Color online) CASE 1: Value at Risk sensitivities. 



of the PDF, reconstructed by Fourier inversion of equa- 
tion (10 1, is fitted by its asymptotic approximation ([6]). 



3.3 The lowest eigenvalue is zero 

The case A* = is discussed here for the set of parameters: 

^ = 0, ^i=l,...,l5 = 1, A* = Ai=i_...^5 = 0, Ai=6, - ,9 = 1 

and A,-io,... ,15 = 2. Figure 5|7 show the same curves de- 
scribed above, again exhibiting full statistical agreement 
with the Monte Carlo simulation. Previous considerations 
about the statistical errors still apply and the only signifi- 
cant difference is the shape of the PDF in Figure [Sj show- 
ing a Gaussian decay of the left tail, as expected from the 
asymptotic expression ([t]). 




0.1 0.5 1 



Fig. 3. (Color online) CASE 1: Expected Shortfall sensitivi- 
ties. 
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Fourier 
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-15 
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Fig. 4. (Color online) CASE 1: left tail of the PDF fitted by its 
analytical approximation corresponding to exponential decay. 




0.1 0.5 



Fig. 5. (Color onhne) CASE 2 (A* = 0): Value at Risk and 
Expected Shortfall. 



3.4 The lowest eigenvalue is positive 

The last case we consider corresponds to the choices: = 
0, (5i=i,...,i5 = 1, A* = Ai^i,...^ = 1 and Xi^5,- ,15 = 2 
(CASE 3). Since A* > 0, the PDF of V is limited from the 



left, see Figure 12 its support is [T4nf,+oo) with T4„t 
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Fig. 6. (Color online) CASE 2: Value at Risk sensitivities. 
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Fig. 7. (Color online) CASE 2; Expected Shortfall sensitivi- 
ties. 
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Fig. 8. (Color online) CASE 2: For A* = the left tail oip{V) 
exhibits a Gaussian decay as expected from ([7|. 



9 — X]fc=i 2A^ ^ —4.75, and the left tail approaches 

with power law scaling {V — Mnf)^^^~^, see equation ([s]). 
As a consequence of the truncation near the origin of the 
V axis, A* and E* are negative for a quite large interval 
oiV* of interest, as illustrated in Figure |9] 
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Fig. 9. (Color online) CASE 3 (A* > 0): Value at Risk and 
Expected Shortfall. The truncation of the PDF near the origin 
is responsible for the negative values for A* and E* over a wide 
interval of V* , see the main text. 
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Fig. 10. (Color online) CASE 3: Value at Risk sensitivities. 



For this set of parameters the unique {5, A) pairs to be 
considered are only {5i, Ai) and {5i, A5), corresponding to 
the sensitivities reported in Figures [lOl ITT] 



4 Conclusions 

In this paper, we addressed the problem of the evalu- 
ation of market risk exposure for non linear portfolios, 
when price changes are approximated according to the 
Delta Gamma Normal model. Exploiting the generalized 
Fourier representation of the probability density function 
for the model, semi-closed form expressions for the Value 
at Risk and the Expected Shortfall were derived. Similarly 
analogous expressions in terms of a single Fourier integral 
were obtained for their first order derivatives; thus our ap- 
proach turns out to be especially useful, since these quan- 
tities are of major practical relevance for the purposes of 
asset allocations and market risk hedging, giving crucial 
information about the sensitivity of the portfolio with re- 
spect to its risky components. It is worth mentioning that 
our approach could be readily extended to compute higher 
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Fig. 11. (Color online) CASE 3; Expected Shortfall sensitivi- 
ties. 
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Fig. 12. (Color online) CASE 3: The PDF is limited from the 
left, decaying as a power law of exponent N/2 — 1 — 6.5. 



order derivatives. All of our formulae lend themselves to 
efRcient numerical evaluation; they have been tested by 
simulating synthetic portfolio scenarios via Monte Carlo. 
This comparison, while confirming the reliability of our 
expressions with a full statistical agreement, highlights 
their potential in practical applications, since full Monte 
Carlo evaluation would require very large samples to ob- 
tain accurate risk estimates. So, the availability of (semi) 
analytical techniques might be welcome. 

A different, yet related, problem stems from dealing 
with time series of limited depth, as it is the case in fi- 
nancial practice. As far as the historical approach is con- 
cerned, this implies a great statistical uncertainty on em- 
pirical quantiles, while for our approach it would be re- 
flected in a noisy estimation of the covariance matrix. We 
did not address this problem in the present work, leaving 
it as a topic of possible future research. A further per- 
spective would be to extend our formalism to different 
financial contexts, where the characteristic function may 
be known but explicit forms for the density function are 
not available, which is especially true for many credit risk 
models. 
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